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We put forward an experimentally feasible scheme for heralded entanglement generation between 
two distant macroscopic mechanical resonators. The protocol exploits a hybrid quantum device, a 
qubit interacting with a mechanical resonator as well as a cavity mode, for each party. The cavity 
modes interfere on a beam-splitter followed by suitable heralding detections which post-select a 
hybrid entangled state with success probability 1/2. Subsequently, by local measurements on the 
qubits a mechanical entangled coherent state can be achieved. The mechanical entanglement can 
be further verified via monitoring the entanglement of the qubit pair. The setup is envisioned as 
a test bench for sensing gravitational effects on the quantum dynamics of gravitationally coupled 
massive objects. As a concrete example, we illustrate the implementation of our protocol using the 
current circuit QED architectures. 
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I. INTRODUCTION 

Experimental creations of nonclassical states of macro¬ 
scopic continuous variable systems have been motivated 
from different perspectives, long-term motivations are to 
explore the limits of the standard quantum mechanics 
and perhaps observe possible corrections at scales where 
collapse phenomena m or gravitational effects become 
significant to fully account for the quantum dynamics in 
table-top quantum optics experiments [5J • A particularly 
interesting case of nonclassical states are entangled states 
between different modes [6]. There are several proposals 
for entangling two mechanical oscillators either in a sys¬ 
tem interacting with a common field CHS or two distant 
resonators without direct interaction mm- Recently, 
a scheme for generating Einstein-Podolsky-Rosen entan¬ 
gled state of two mechanical systems has been proposed 
m- These schemes mainly deal with creation of bipar¬ 
tite Gaussian entanglement demonstrated by the second 
moments of the phase space quadratures. 

Heralded entanglement generation is tailored to en¬ 
tangle two remote particles where the entanglement is 
produced conditionally based on measurement outcomes 
without any direct interactions between the particles m- 
This method has been proposed for entangling distant 
qubits and recently has been implemented in a num¬ 
ber of different systems cansi. It would be intrigu¬ 
ing to extend such technique for creating remote entan¬ 
glement between macroscopic mechanical resonators. In 
this work, we propose a protocol for generating entan¬ 
gled states of two well-separated noninteracting macro¬ 
scopic mechanical resonators. Our heralding technique 
is useful for preparing entangled coherent states m , 
i.e. A/±(|a)|/3) ± |/?)|a)), between massive objects (see, 
Fig. 0 indicating entanglement between first moments 
(centers of mass) of the resonators. This type of contin¬ 
uous variable entanglement is non-Gaussian which is— 
unlike Gaussian entangled states—characterized by neg¬ 
ative Wigner function in phase space representation, and 


thus is fundamentally inconsistent with any classical de¬ 
scription. An interesting scenario here could be to con¬ 
sider the entanglement between two test masses which 
are only gravitationally interacting with each other and 
to track a genuine gravitational decoherence on the en¬ 
tanglement dynamics |20| . Furthermore, the setup is use¬ 
ful for Bell test performed on remote macroscopic sub¬ 
systems (3EDE31. We show that the proposed protocol 
is experimentally feasible with the current technology of 
the so-called circuit QED devices [23] . 


II. MODEL 

We use a particular hybrid quantum device for imple¬ 
menting our entangling protocol. The hybrid device is 
composed of a mechanical oscillator with frequency w m , 
a photonic cavity mode with frequency u> c , and a qubit 
system made up of ground | g) and excited |e) energy lev¬ 
els. We also need an extra, auxiliary, excited state |/) for 
conditionally entangling the qubit with an emitted cav¬ 
ity photon. Solid-state qubits such as spin of a nitrogen- 
vacancy center p?5] or superconducting transmon qubits 
[261 have such a structure. The qubit strongly interacts 
with both the mechanical and the electromagnetic modes, 
i.e., its coupling rates to these modes are greater than the 
respective decoherence and damping rates of the system. 
The mechanical resonator is pre-cooled to its motional 
ground state which can be achieved by sideband cooling 
of the mechanical resonators via its coupling to the qubit 
m- Since typically the thermal excitation numbers of 
the cavity and qubit itself are very small at cryogenic am¬ 
bient temperatures, this way one could practically pre¬ 
pare the system very close to its ground state. In order 
to create a mechanical entangled coherent state we are 
specifically interested in a qubit-mechanics interaction of 
the form Xjx\j){j\ with j = {g, e, /} and x being po¬ 
sition of the resonator. This describes a state-dependent 
force generating a time evolution of state-dependent dis- 
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distinguishability of the two mechanical coherent states. 
Thus, one needs to increase it by enhancing the qubit- 
mechanical coupling rate which is actually very demand¬ 
ing from a technical point of view though notable efforts 
for attaining it [25] ■ Instead, the long coherence time of 
the current technology of qubits allows for exploiting the 
method proposed in Ref. [33] to increase this distance. 
The idea is to apply a sequence of [Tr] g ++ e pulses to flip 
the qubit periodically synchronized with the resonator 
frequency. By choosing an odd number of such pulses 
at half mechanical period time intervals it/tjJ m the state 
prepared at the end of this stage will be -jA\g,a,Q>) + 
|e, -a,0)]A(s), where a = (N p + l)(A e - A g )/w m with N p 
the number of pulses. 

(3) Conditional photon emission. We now ‘copy’ the 
qubit excitations into the cavity photons by employing 
their third state |/). One first applies a [ 7 r]^/ pulse 
to bring up the qubit from |e) to |/). Once being back 
in |e), due to relaxation or other mechanism, a single 
photon will be emitted conditioned on the state of the 
qubit prior to the [jr] e ++f pulse. This, therefore, results 
in the following state at each party 


FIG. 1. (color online) (a) and (b) Schematic representation 
of different types of entangled cat states of two suspended 
masses moving harmonically along rr-axis. a is the displace¬ 
ment amplitude from the equilibrium position, (c) Sketch of 
a possible experimental implementation with quantum elec¬ 
tromechanical circuit. 


placement of the resonator. In this model the cavity 
mode resonantly interacts with the qubit described by 
the Jaynes-Cummings Hamiltonian \(\f){e\a + \e)(f\al), 
where a is the cavity field annihilation operator. This 
leads to oscillatory |e) O |/) transition which can be 
exploited to entangle the qubit with a traveling photon. 


III. PROTOCOL 

Having the above system in two remote sites, the fol¬ 
lowing protocol can be applied for entangling their dis¬ 
tant mechanical resonators: 

(1) Initialization. The separated parties A and B 
prepare their system in the ground state: \g,0, 0)a <§) 
|3;0i0)b) where within each ‘kef the first argument 
stands for the qubit, the second for the mechanical mode 
being in ground state, and the third argument denotes 
the vacuum mode of the cavity. Immediately after¬ 
wards a [7r/2] g< _ ) . e pulse is applied to each qubit giving 
|+,0,0)a® |+,0,0) b where |±) = -A=(\g) ± |e)). 

(2) Displacing the resonators. After the state ini¬ 
tialization, the qubit and its respective mechanical mode 
interact for a time duration of r. For a coupling rate 
A j the maximum achievable displacement is ay = — T 2 - 
which is obtained when r = 7r/w m . Therefore, the me¬ 
chanical modes will be conditionally displaced in their 
phase space. The distance \a e — a g \ determines the 
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[|.g, a, 0) -(- |e, 


a ’ 7 ] A(B) ’ 


(1) 


already indicating a local three-body entangled state be¬ 
tween, the qubit, mechanical resonator and cavity pho¬ 
ton. 

(4) Heralded hybrid entanglement. The cavity 
modes of the parties A and B interfere on a 50:50 beam¬ 
splitter resulting in 
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|.9, a)A|e, ~a) B + |e, -a) A \g, a) B 
\g,a) A \e,-a) B ~ \e,-a) A \g,a) B 
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|1,0)d 
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where | , ) B = | )dx \ )d 2 denotes the state of the output 
modes after the action of the beam-splitter. An essen¬ 
tial feature appearing here is Hong-Ou-Mandel effect in 
the second line. This effect has been recently realized 
in microwave regime m This arises due to the indistin- 
guishability of the two input mode photons interfering on 
the beam-splitter. Therefore, generation of indistinguish¬ 
able photons at the input ports is crucial for successful 
heralded entanglement. The first two lines show no en¬ 
tanglement while the third and forth lines contain hybrid 
entanglements. By inspecting the above expression, we 
realize that the projection onto an entangled state in a 
single shot measurement can be achieved by adopting a 
suitable photon detection scheme. Here, we choose it 
to be the photon-number parity detection at the beam¬ 
splitter output modes. This can be done by placing a ‘de- 
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tector’ qubit at each output port j £ {1,2} and perform¬ 
ing a controlled-phase gate C n = \g){g\ + e 1 ^ a o j aD j | e ^ e | 
causing a bit flip |±) —> |=f) in the corresponding qubit 
conditioned on the arrival of odd number (here, only one 
photon) of photons. Therefore, we can map the photon 
parity of the output modes onto the qubit states and 
serve as parity detectors. The bit flip can be detected in 
a Ramsey measurement by inserting C n between two ir/2 
pulses on the detector qubits. Accordingly, detection of 
a bit flip in either output port heralds the corresponding 
qubit-mechanical entangled state 

\g, cx)A\e, ~u)b ± je, -a) A \g,a) B 

^2 U 

The above state involves superposition of the two res¬ 
onators being in different relative distances (see, Fig. &)• 
Another form of entangled state can be created by sim¬ 
ply overturning the direction of the state-dependent force 
in one of the qubit devices giving superposition of the 
center-of-mass in different locations (see, Fig. §>)• 

(5) Entangled cat state. Finally, in order to disentan¬ 
gle the mechanics from the qubits and obtain a purely 
mechanical entanglement, each site applies a [7r/2 ] g <^ e 
pulse to its qubit resulting in: 

{\9)a\9)b ~ \e) A \e) B )W+) + (|e)^|g)B ~ Iff) A\e) B )\j>-) 

2\[2 

Here, we have defined the state \ip±) = \ — ol)a\oi)b ± 
\ol)a\~&)b f° r the mechanical parties which already ex¬ 
hibits entanglement between the first moments of the 
two distant mechanical resonators. The users now can 
read out their qubit and post-select \i/j±) with probabil¬ 
ity p± = |(1 ± ) according to the measurement 

outcomes. 


IV. IMPLEMENTATION 

Here, we focus on a specific implementation of our 
protocol. High controllability and tunability of super¬ 
conducting qubits makes them a versatile tool for engi¬ 
neering different regimes of interactions and control over 
photonic and mechanical systems |i 31I 132) . There are im¬ 
portant achievements in fabricating such hybrid devices 
where a vibrational mode of a mechanical resonator is 
coupled to a superconducting qubit and at the same time 
the qubit strongly interacts with a coplanar microwave 
resonator [25} [35]. These make the hybrid circuit quan¬ 
tum electrodynamical devices a promising framework for 
implementing our protocol. We consider, in particular, a 
recently proposed hybrid electromechanical circuit which 
can basically be fabricated and employed by the current 
technology [53] . 

The device consists of a transmon qubit capacitively 
coupled to a microwave coplanar waveguide and a me¬ 
chanical resonator (see Fig. [Tj) . The reduced anharmonic- 
ity of a transmon qubit makes it possible to access its 


higher levels, thus giving us the three-level ladder system 
for producing cavity photons from the qubit excitations. 
The Hamiltonian of the system at site A (similarly at site 
B) reads [33;, 331 

H s = u> c a\aA + * v m (pA + x 2 A ) - E 3 cos <^4 

+4E c (x A )[nA - v(aA + a^)] 2 , (3) 

where, Ej and Ec are Josephson and charging energy 
of the superconducting qubit, respectively. Here, < 2,4 is 
the annihilation operator of the microwave photons inside 
cavity A, while xa and pa are respectively the normal¬ 
ized mechanical position and momentum operators with 
commutation relation [xa,Pa] = *• Also, tia and ip a are 
the superconducting charge number and phase operators, 
satisfying the commutation relation [ip a , ua] = i and rj 
is the Lamb-Dicke parameter. Since the charging en¬ 
ergy of the qubit depends on the position of the mechan¬ 
ical resonator, one arrives at A^^-vYx for the transmon- 
mechanical mode interaction featuring state-dependent 
force on the mechanical mode. Here, fi = Yin n\n)(n\ 
where n is the number of exchanged Cooper pairs with 
corresponding eigenstate |n). By writing this in energy 
eigenstates of the transmon Hamiltonian one retrieves 
the interaction given above. The influence of the num¬ 
ber of transferred cooper pairs by the cavity field also 
results in a Jaynes-Cummings interaction between the 
qubit transitions and the cavity mode with Rabi fre¬ 
quency x = 8Ecrj(f\n\e} [see Appendix A], This can be 
employed for coherently converting the qubit excitations 
into the cavity photons as explained below. 

In such devices, an external magnetic field applied to 
the Cooper pair box tunes the transition frequencies of 
the qubit. Thus, bringing either of the transmon tran¬ 
sitions into resonance or taking them off-resonance from 
the cavity mode frequency. In the third step of the pro¬ 
tocol, we specifically are interested in the situation where 
frequency of the cavity matches qubit’s first to second ex¬ 
cited state transition: w c = H/ —fl e . To make a travelling 
microwave photon conditioned on the state of the qubit, 
one first applies a 7r-pulse which flips the qubit from |e) 
to I/) then takes the state transfer interaction between 
transmon and cavity into resonance for half period of a 
Rabi oscillation 7 r/ 2 y. This brings the qubit into its first 
excited state accompanied with emission of a microwave 
photon. We remind that high fidelity single-qubit gate 
operations can be performed by properly shaped disper¬ 
sive microwave pulses which can be used for controlled 
flipping of the qubit in either of its transitions [35} [37]. 

Emitting indistinguishable photons from each cavity is 
essential for a faithful projection onto an entangled state. 
Indistinguishability of the photons can indeed be guaran¬ 
teed by employing two low finesse microwave resonators 
on the sites. This leads to a broad wave-packet which 
increases overlap of the photonic wavepackets, and there¬ 
fore, their indistinguishability. Its other consequence is 
lowering emission time of the photons which basically 
makes them travelling photons. Therefore, reducing the 
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protocol run time. 

Moreover, the controlled phase gate performed on de¬ 
tector qubits discussed in the fourth step of the protocol 
can also be realized using superconducting qubits dis- 
persively coupled to a coplanar microwave resonator at 
the output of the beam-splitter. The Hamiltonian of this 
part of the system is that of a qubit dispersively coupled 
to a cavity with frequency ui c , the same frequency of the 
site cavities (see below) 

2 2 

H p = (w c + e][e\)a ] D .a Dj + (H e + ^)|e)(e|, (4) 

where A = fl e — w c is the qubit-cavity detuning. Because 
of this dispersive coupling the unitary time evolution op¬ 
erator of the detector qubits in a frame rotating at the 

x 2 

qubit and cavity frequencies is exp{z -£* l e )( e l}- 

Our goal is to have a conditional 7r-phase shift given that 
the photon number in the cavity is odd. Therefore, the 
interaction time of the qubit and the photons must satisfy 
Xpt/A = ir. By storing the incident photons in a high 
finesse cavity the interaction time between photons and 
the detector qubits can be increased sufficiently large. 
Therefore, after arrival of the photons to these secondary 
cavities one waits for nA/Xp seconds then measures the 
parity by performing a Ramsey pulse sequence on detec¬ 
tor qubits [38] . Today’s experiments are able to measure 
the parity of the storage microwave cavity with fideli¬ 
ties above 90% via another read out cavity coupled to a 
transmon superconducting qubit [231 . 

It worth mentioning here that to maximize the absorp¬ 
tion of these photons (after mixing at the beam-splitter) 
into the detecting cavities two conditions must be met: 
First, the resonance frequency of the cavities must match, 
their linewidth must be close to each other. The former 
condition, in principle, is easily met by employing equal 
frequency cavities, while the latter sounds contradictory. 
On the one hand, in the detection parts, we need to em¬ 
ploy high finesse cavities in order to store the arrived 
photons and give them enough time to rotate the qubit 
states. On the other hand, as discussed above, site cavi¬ 
ties must have high decay rates to ensure indistinguisha- 
bility of the outgoing photons. This can be resolved by 
employing a cavity with tunable decay rate in the detect¬ 
ing parts [32] ■ These cavities are designed such that can 
be tuned in situ to different decay rates for the purpose 
of maximal capture, storing, and retrieving microwave 
photons. In our case, the detecting cavities can be first 
tuned to a linewidth which allows for maximal capture of 
the incoming photon(s), then the photon(s) can be stored 
for performing the parity measurement. 


V. DISCUSSION 

The question of how to verify the nonlocal coherence 
involving the mechanical parts can be approached by 
looking into the dynamics of the two qubits coupled to 







FIG. 2. (color online) (a) Periodic revival of the concurrence, 
undergoing decoherence, signifies the entangled superposition 
between the masses for small (dashed line) and large (full line) 
displacement amplitudes. For larger amplitudes the decay of 
entanglement revival under decoherence occurs with a higher 
rate, (b) Preparation of large entangled superposition using 
[7r] 9< ^ e pulse sequence, and monitor the concurrence after¬ 
wards without decoherence (full line) and with decoherence 
(dashed line). 


their respective resonators. Since ±a in Eq. ([2]) is pe¬ 
riodic in time it becomes zero at some points at which 
the mechanical part is disentangled from the qubit part, 
and thus leaving the qubit pair in a maximally entangled 
state. This is manifested in the entanglement revival in 
the qubit pair which can be considered as a signature 
of the entangled superposition state between the masses 
and qubits, i.e., Eq. ([2]), at intermediate times. A sim¬ 
ilar consideration has been used previously for probing 
macroscopic superposition states in Refs. [4H1I41] . In our 
case this feature can point to even a stronger indication 
of nonclassicality. Monotonicity of the entanglement al¬ 
lows the entanglement revival if and only if the qubit 
pair has access to a global coherent operation (e.g. di¬ 
rect interaction) inducing this entanglement. Since there 
is no direct interaction between the qubits this revival 
has to be provided by the dynamics of the resonators be¬ 
ing in an entangled superposition with the pair prior to 
the revival time. This is an unambiguous way of veri¬ 
fying entangled superposition including the mechanical 
parts. Because the action of separable resonators on the 
qubits fall under the local operation and classical commu¬ 
nication (LOCC) operations which cannot increase the 
entanglement. For this purpose, one should use an en¬ 
tanglement monotone such as concurrence [42] or nega¬ 
tivity [43] for monitoring the entanglement. Concurrence 
has already been measured experimentally on supercon¬ 
ducting qubits with high fidelity using state tomography 
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[331 [IS]. In Fig. [ 2 ] we plot time evolution of the con¬ 
currence of the two qubits for two different displacement 
amplitudes and the effect of decoherence on this evilu- 
tion. The detailed analytical and numerical analyses are 
presented in Appendix B. 

Information transmission about mass’ positions due to, 
for instance, the interaction with the environment sup¬ 
presses the magnitude of the entanglement revival and 
the state evolves towards statistical mixture with the 
same rate regardless of which form of mechanical entan¬ 
glement (a) or (b) shown in figure [l] was created [see 
Appendix C]. However, if we think of an unconventional 
situation in which the two masses gravitationally inter¬ 
act with each other, then the masses’ configuration in 
the mechanical parts of the entanglement becomes im¬ 
portant. One model of gravitational decoherence is pre¬ 
sented in Ref. [35] which is equivalent to a different for¬ 
mulation proposed in Ref. [20]. In these models the deco¬ 
herence rate is completely determined by the gradient of 
the gravitational force between the masses. The gradient 
of the gravitational forces in the two types of entangle¬ 
ment shown in Fig.Qa) and (b) are different. Therefore, 
adding this new source of noise may lead to a detectable 
gap in the decay of entanglement between the two cases 
undergoing identical environmental decoherence, yet dif¬ 
ferent gravitational decoherence rates. This gap then can 
be attributed to a genuine gravitational effect. This is an 
intriguing strategy. Because, the main challenge in ob¬ 
servation of gravitational effects is that it is hard to dis¬ 
tinguish the intrinsic gravitational field effects from those 
of environmental (conventional) decoherence in quantum 
dynamics as both these sources of noise leads to similar 
reduction of the superposition states. 


VI. CONCLUSION 


and M.A. by the Alexander von Humboldt Foundation 
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Appendix A: Implementation Hamiltonian 


The Hamiltonian describing the device we are consid¬ 
ering in this paper is composed of the transmon qubit 
coupled to the mechanical resonator via its charging en¬ 
ergy, that is (H = 1 ) 

H s = uJcC^a + ^ ui m (p 2 + x 2 ) — Ej cos(7r<lJ/<ho) cos<£ 
+4Ec(i) [n — n g — rj(a + a 1 ’)] 2 , (Al) 


C 1 

where 77 = -g i ( Huj c /2C )$ is the circuit QED Lamb-Dicke 
parameter and n g is the induced dc gate charge. Here, >I> 
is the externally applied magnetic flux through the super¬ 
conducting loop and $0 = h/2e is the superconducting 
flux quantum. The operators h and (p denote the number 
of Cooper pairs transferred between the islands and the 
gauge-invariant phase difference between the supercon¬ 
ductors, respectively. 

Since in transmon qubits Ej/Eq 1 , the nonlinearity 
of the qubit is reduced such that the higher levels play a 
role in its dynamics. We simplify the above Hamiltonian 
by first Taylor expanding the Eq (£) and keeping only to 
the first order in x. Then applying rotating wave approx¬ 
imations and truncating the transmon Hilbert space to 
its first three levels. This brings us to 


H s = w c a f a + -w m (p 2 + x 2 ) + + Xjx)\j)(j\ 


+X 


LV^kXsl + |/)(e|)a + H.c. 


(A2) 


In this work we proposed an experimentally feasible 
protocol for creating heralded entangled cat states be¬ 
tween spatially separated mechanical resonators. The 
scheme can be implemented in currently available cir¬ 
cuit QED architectures. We expect that the setup will 
provide a useful platform for experimentally probing the 
interface between gravity and quantum physics where the 
mechanical resonators in the present scenario serve as 
test masses undergoing gravitational decoherence. Fi¬ 
nally, we think that going beyond a single particle super¬ 
position, and monitoring the dynamics of ‘gravitation¬ 
ally’ different types of entangled states may provide a 
more visible test of gravitational decoherence models in 
a new regime of gravitational sensing. 
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with j = {<?,e, /}. Here, a is the annihilation oper¬ 
ator of the microwave photons inside the cavity, while 
x and p are respectively the normalized mechanical po¬ 
sition and momentum operators with commutation re¬ 
lation \x,p\ = i. The Hamiltonian (A2|, already fea¬ 
tures a state-dependent force on the mechanics via the 
qubit which is crucial in our protocol. The second line 
of the Hamiltonian is the generalized Jaynes-Cummings 
transmon-cavity interaction. 

The probe qubits discussed in the fourth step of the 
protocol could be realized by superconducting qubits dis- 
persively coupled to the coplanar microwave transmission 
line at the output of the beam-splitter. Here, the goal is 
to have a parity flip in the qubits conditioned on the odd 
incident photon numbers. Therefore, the interaction time 
of the qubit and the photons must satisfy y p f/A = 7 r. 
The pulse duration of the incident single-photon states 
is roughly the same as the time it has taken to be emit¬ 
ted in the sites, i.e. 7 r/ 2 y. Therefore, one needs to fulfil 
Xp = 2yA for getting a half rotation about the 2 -axis for 
every photon. Since we have A y p from the dispersive 
coupling regime, this can be achieved only for y p y. 
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This looks, however, impractical for a waveguide because 
of the finite coherence time of the qubits. By storing 
the incident photons in a high finesse superconducting 
resonator the interaction time between photons and the 
probe qubits can be significantly increased. 


Appendix B: Monitoring the entanglement dynamics 
under the influence of decoherence 


describes qubit dephasing acting locally with dephasing 
time T 2 = l/ 7 q and 


= ^{n + 1 )( 2 bjpb] - b]bjp - pb\bj) 

+ yli(2 b)pbj - bjbjp - pb]), (B3) 


We aim to monitor the entanglement dynamics of the 
two-qubit system in which the qubits and their respec¬ 
tive resonator evolve under the influence of the interac¬ 
tion with the environment. Here, we use the Woofer’s 
concurrence, which is defined as C{p) = max{0,v / W — 
— VX 3 — vXt} where Xi are the eigenvalues of the 
matrix pp = pa y <S> a y p*a y (8 <J V and Ai is the maximum 
eignevalue. 


describes the mechanical dissipation within a single res¬ 
onator, where = oj m /Q m is the mechanical damp¬ 
ing rate for mechanical resonator with quality factor Q m 
and n = 1 / [e huim / k ^ T — 1 ) is the equilibrium occupa¬ 
tion number (identical resonators and dissipation is as¬ 
sumed). In the high temperature limit ii> 1 we obtain 
T th = q m n ~ k^T/(hQ m ) as the relevant mechanical 
decoherence rate. 


1. Environmental decoherence: Analytic treatment 

As an illustration, in this section we treat a simple sce¬ 
nario of probing the entanglement dynamics under deco¬ 
herence in which the qubits coupled with constant cou¬ 
pling strengths to their respective mechanical resonators. 
We evaluate the effect of qubit decoherence and mechan¬ 
ical decoherence due to a weak coupling of the resonator 
to a finite temperature bath where the Markovian mas¬ 
ter equation of a Lindblad form is applied. Therefore, we 
model these decoherence processes by a master equation 
of the form 


Let us now consider the effect of decoherence on the 
entanglement dynamics of the qubit pair after the prepa¬ 
ration of the hybrid entangled state p{t 0 ) obtained by 
the heralded detection. The state at later time t after a 
period of r = t — to obtained from Liouvillc superoper¬ 
ator of the time evolution p(t) = e CT p(to) which is the 
solution to Eq. (Bll. In particular we are interested in 
the dynamics of the off-diagonal term of the qubit pair’s 
reduced state 


Pge,eg(t) = {g\{e\p{t)\e)\g) = ( g,e\p qAqB {t)\e,g)p^'2m B {t)■ 


P = C-P = ~i[H,p\ + Y £qjP+ Y Cm i p ' ( B1 ) 

j=A,B j=A,B 


where 


C qjP = W ( Y \ k )( k \jP\ k )( k \j ~ P 


(B2) 


Now, let us analyze a scenario in which after a fast prepa¬ 
ration of the hybrid entangled state at time t 0 the system 
undergos decoherence processes according to Eq.(Bl| at 


the later time t (t to) For the expectation value at 
t = to + t we need to solve 


Pge,eg(t ) = ( 5 , e\ (e Cr \g, e)(e, g\(g, e\p qAqB (t 0 ) |e, g)p 9 ^?a B (*o)) |e, g) ■ (B4) 

According to the master equation ( |B1[ ) and Hamiltonian (1) in the main text this operator evolves as 

Pge,eg — ^ ~ [bjbj , Pge,eg\ T £"mj Pge,eg T ^Qj Pge,eg J* 2 ^e,A (b^ T ^a) Pge,eg T ^A e,B{bj^ T ^s) Pge,eg• (B5) 

j=A,B 

We have assumed \ g = 0 for both of the systems. 

From above we can obtain the time evolution of an entanglement monotone for the two-qubit system. In the special 
case we consider, concurrence simplifies to 

C(p{t)) = 2 |T rm A ,m B {Pge,eg{t)} \ ■ (B6) 

Note that, the concurrence is completely characterized by the off-diagonal term of the qubit pair’s state in which the 
qubits undergo only dephasing and the mechanical resonators interact with a finite temperature bath. The plot is 
shown in Fig. 2(a). In the absence of the resonator decoherence 

Tr m A ,m B {Pge,eg(t)} = ^' ,T (^A [ct e (7")]) mA (Ps [a e (r)]) mB . 
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The displacement amplitude is a e (r) = X e /uj lrl {e~' lulrnT — 1) which is periodic in time, and thus the concurrence 
demonstrates collapse and revival. Here, we have taken A eyA = ^e,B = A e . Concurrence is the same regardless of 
which state 


or 


1 

71 


i 

71 


(I g,cx) A \e,-a) B ± \e,~a) A \g,a) B ) 


( \g, -a) A |e, -a) B ± |e, a) A \g, a) B ) 


(B7) 


(B8) 


was prepared at time to- For entangled state (B 8 ) the exerted force from the second qubit is reverted with respect 
to the x-axis. Therefore, the associate configuration indicate superposition between two different relative distance 
between the resonators. 

t) = 


Equivalently, we can define the Wigner characteristic function xTg A {P-A,t) = (U(f3 A ))m A and 
(V(/3 B )) mB and write 




Tr m A m B {Pge,eg{t)} — zt ^ q Xge(PA:t)Xeg{pB-,t). (B9) 

The evolution of the characteristic function is given by the Fokker-Planck equation 

Xeg(P) = i - A*/?* XegW) ~ ^ (2« + 1) \P\\eM (BIO) 

+ (^f~) XeflO 3 ) - iX « (jp ~ Xeg(P): 

where B = u> m + * 7 m / 2 . We solve this equation in three steps. First, we make the ansatz 

Xegtf, t ) = e ^W e ^*W-^«W Xl (/3, t), (Bll) 


where k = —?Bk — iX e /2, and <f> = —A e («(t) + For the remaining equation for X j(/3,t) we have 


d 


d 


Xi(P,t) = M (0/3- A e ) — - (B*/3* - A e ) — X/ (/3,i) - -^-(2n + 1)|/3| X /(/3, t). 


a/3 




We now make the second ansatz 


X i(/S,t) = e -(«+i)(l/3| 2 -^*(t)-rs( t )+C W) X n(p,t), 


(B12) 


(B13) 


where k = —/Bit — *A e , and £ = —i\ e (k — k*). This leaves us with the remaining equation for Xn{i 3, i)> which is given 
by 


Xn{P,t) = 


*(0/3- X e )-~i(n*(3*-X e ) — 


a/3 


a/3* 


xn(P,t)- 


(B14) 


This equation is solved by a function of the form 

A.7/03,i) = X jj fa = e lQt p - iX e e lQs ds^J , (B15) 

and the specific expression for Xii( x ) is determined by 
the initial conditions 

Xii(x) = e^ m+1 ^ 2 / 2 Xeg {x, t = 0). (B16) 

For an initial thermal state Xii( x ) = 1 and therefore 

Xfl eG8 = 0,t)=e-( fi +^«t). (B17) 


In the limit 7 m <C oj where the dissipative part of the 
master equation is valid we obtain 


at) ~ 


2a: 


2 r 


(1 - cos (u m t)e 7mt/2 ) + ^7 


(BIB) 


This shows that for A e ~ the signal of a sin¬ 
gle measurement decays with a total decoherence rate 
Tdec = (2n +1) 7 m - For k B T huj m we obtain n 7 m ~ 
k B T/HQ m . Similar conclusions are obtain, when starting 
from a precooled state or for A e <C ui m , when a 7 r-pulse 
sequence is obtained to amplify the displacement ampli- 









tude. 

We should remark that in the decoherence model it 
is assumed that the oscillator damping is very small 
( 7 m ujm)- We use this model of damping in this section 
to provide a simple illustrative evaluation of the effect of 
the resonator’s dissipation on the qubit pair’s entangle¬ 
ment dynamics. In the next section we use quantum 
Brownian motion master equation giving more accurate 
results in higher damping rate and reduce to above mas¬ 
ter equation using rotating wave approximation. In the 
next section, we treat this scenario numerically in a more 
involved decoherence processes. 


2. Numerical simulation of the full dynamics of the 
system 


The system dynamics is composed of free evolution and 
the dissipations. The free evolution part is captured by 
the system Hamiltonians in the mechanical sites given by 
© and (3) in the main text, and the dispersive Hamil¬ 
tonian describing the dynamics of the probe sites given 
by the Eq. (4) in the main text. And the dissipations and 
decoherence stemming from dephasing of the transmon 
qubits and relaxations in the cavity modes, and finally 
thermalization in the mechanical modes. Here, we as¬ 
sume that the cavity decay rates is majorly due to its 
coupling the the transmission lines. This will not lead 
to serious restriction as it is almost the practical case. 
Moreover, note that the cavities are irreversibly coupled 
to each other via transmission lines mediated by a beam¬ 
splitter. These are the major sources of imperfection and 
one includes them all in a single master equation to study 
full dynamics of the whole system. 

P = —i[H, p] + (£<j + £-m + L c )p, (B19) 


where H = H s< a + H S)B + Hp.o, + -ff p ,_D 2 is the total 
Hamiltonian which is composed of the both mechanical 
sites and the detection facilities denoted by ‘p’ subscript. 
We have also decomposed the dissipator into three part 
each of which for a distinct form of energy: photonic, 
charge, and phonoic. The mechanical resonators damp¬ 
ing rate is 7 m and because of the low frequency its dif¬ 
fusion is largely affected by the environmental thermal 
phonons n. The mechanical dissipator is already intro¬ 
duced in Eq. (|B3|). 


For the qubits, the dissipator is can be divided into the 
four local dissipators (for four qubits) C q = La + L B + 
Ljj 1 + Ljj 2 . This assumption is true so long as there is no 
dissipative coupling between the qubits, and keeps to be 
true in our case, since the qubits are well separated from 
each other and any potential coupling between them will 
happen coherently via the coupled cavities (see below). 
The dissipator of every qubit must include both relax¬ 
ation and its pure dephasing. For example La is given 


by 

La[p] = ( 2 \j)( k \p\ k )U\ ~ \ k )( k \P ~ P\ k )( k \) 

j<k 

+^E(b'XilPb') 01 -p)- (B 20 ) 

3 

where |fc) with k = {g, e, /} are the three lowest trans¬ 
mon states. The first of the above dissipator and the 
second line correspond respectively to the relaxation and 
the pure dephasing of qubit A. The relaxation happens 
with rate 7 a and the pure dephasing rate is 7 a such that 
the total dephasing time of the qubit is 1/T| = ja +7a/2. 
The same arguments hold for the remaining three qubits. 

Finally, for the cavities, the photons experience both 
decay into the coupled transmission line and loss (ab¬ 
sorption and decay into the free space), leading to the 
total decay rate K to t = K; n + Ki oss . However, the cur¬ 
rent technology superconducting resonators have negligi¬ 
ble loss Ki n ^ ki oss and the microwave photons mostly 
escape to the coupled transmission line Ktot ~ Kin- We 
therefore, have a dissipator in the Lindblad form with 
the decay rates k s and k p for the site and probe cavities, 
respectively. We also include the unidirectional coupling 
of the detecting cavities and the site cavities. Actually, 
this coupling is not direct and it happens via the beam¬ 
splitter. Thus, the detecting cavities are irreversibly fed 
by both of the site cavities. Since we are considering a 
50:50 beam-splitter, the effective mode which would cou¬ 
ple to the L>x and D 2 cavities are (ua ± 0 , 3 )/V 2. By 
taking this into account, the following Liouvillian holds 
for the cavity modes 

£cp = |^-(Ta + L b ) + -L(Lc 1 + Tc 2 )|p 

)p] + . 

+ [4> 2 .(^W] + (B21) 

where Lap = 2 opo 1 — o*op — po^o is the Lindblad dissi¬ 
pator. Here, the last two lines express the irreversible 
coupling between the probe and site cavities including 
the beam-splitter mixing effect and e is the efficiency of 
the transmission channels which takes values 0 < e < 1 
mm- The parameter e contains the waveguide losses, 
which are typically negligible, and reflection of the mi¬ 
crowave photons at the port of the detecting cavities and 
the beam-splitter. As we have discussed in the main text 
the reflection effects can be minimized by appropriate 
choice of devices, therefore, giving e very close to one. 


3. Simulation’s parameters 


One numerically solves the full master equation (B19) 
with the above Liouvillians and by a post-selection sim¬ 
ulated by a projective measurement one of the states 
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to) 



FIG. 3. (a) Time evolution of entanglement of two qubits in terms of concurrence with (dashed line) and without (full line) 
decoherence, and (b) probability of finding the A qubit in its ground state and the B qubit in its firs excited state. This 
quantity is obtained by simulating the protocol and the evolution of the prepared state using the parameters listed in Table [I] 


(2) will be obtained. However, this cannot be done by 
the available computational resources because of the very 
large system size. Instead, we turn to simulate the pro¬ 
tocol step by step. Therefore, we first solve the master 
equation for each local site coupled the waveguide at the 
output which at the end of the third step of the proto¬ 
col gives Qj with j = A,B, then we merge their photon 
parties by assuming a perfect 50:50 beam-splitter. In the 
next step, the dynamics in the probe cavities are simu¬ 
lated with the initial separable quadripartite qubit-cavity 
state £>o<S)|+X+|di®|+X+|d 2 , where g 0 = Tr 9jm { U(g A ® 
g B )U^} is the photonic party of the state mixed at the 
beam-splitter. Here U = exp{ * (d^d/j — a^a^)} ex¬ 
presses the beam-splitters unitary operation, while gA 
and gs are the outputs of the first stage of the simula¬ 
tions. The measurements are simulated as perfect pro¬ 
jections. Therefore, the post-selected state according to 
the parity of the detecting qubits and state of the site 
qubits gives us the final state which turns out to be an 
entangled coherent state. 

In Fig. [3] we plot the evolution of the probability of 
simultaneously finding the A and B qubits in the ground 
state and excited state, respectively. The parameters 
used here are feasible with the current technology [see 
Table [I] . In the detection section, a microwave resonator 
with k p /2tt = 20 kHz is considered. The detector qubits 
are transmon qubits with the same properties as the site 
qubits and operated at Ej/Eq = 40 which makes them 
well away from the cavity resonance (x P /A « 0.02). This 
means having nA/x p ~ 1/4k s giving enough time for 
changing the qubit parity. The ambient temperature is 
taken to be T = 25 mK. 

It worth to mention here that in the detection parts, 
we need to employ high-Q cavities in order to store the 
arrived photons and give them enough time to rotate 
the qubit states and thus perform the parity measure¬ 
ment. On the other hand, the photons leaving the A 
and B cavities must have broader band to ensure indis- 
tinguishability, the key point for Hong-Ou-Mandel effect. 
To maximize the absorption of these photons (after mix¬ 
ing at the beam-splitter) into the detecting cavities two 


TABLE I. Parameters of the system. 


Quantity 

Symbol 

Value 

Mechanical mass 

m 

3 Pg 

Mechanical frequency 

0Om/2n 

1 MHz 

Mechanical quality factor 

Qm 

10 5 

Transmon-mechanics coupling rate 

Xe/2-K 

50 kHz 

Josephson energy 

Ej/2n 

35-55 GHz 

Charging energy 

Ec/2n 

0.5 GHz 

Transmon relaxation rate 

7?/ 2 7r 

5 kHz 

Transmon pure dephasing rate 

7 9 /27t 

20 kHz 

Transmon-cavity coupling rate 

x/2tt 

45 MHz 

Cavity frequency 

uj c /2ty 

11 GHz 

Cavity decay rate 

k s /2iy 

200 kHz 


conditions must be met: First, the resonance frequency 
of the cavities must match, their linewidth must be close 
to each other. The former condition, in principle, is eas¬ 
ily met by employing equal frequency cavities. However, 
the second condition requires n p « k s , which is in con¬ 
tradiction with our requirements of the parity measure¬ 
ment. This can be resolved by employing a cavity with 
tunable decay rate in the detecting parts [39]. These are 
designed such that can be tuned in situ to different decay 
rates with three orders of magnitude difference. In our 
case the difference between k s and « p is only one order 
of magnitude: 1/k s ss 8 x 10“ ' s and 1/k p «8x 10 -6 s 
within the reported cavity lifetimes in Ref. [3TJ|. 


Appendix C: Gravitational decoherence of collective 
modes of two gravitationally coupled harmonic 
oscillators 

For the sake of gaining more insight into the idea dis¬ 
cussed in the main text let us focus on a simple and 
ideal example. An interesting scenario happens if the 
two mechanical resonators are gravitationally coupled to 
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each other for sufficiently large gravitational coupling 
strength. In this case dynamics of the two resonators 
is described by two independent collective modes of os¬ 
cillations with coordinates x + = (xa + xb)/V%, for the 
center of mass mode, and £_ = (xa — Xb)/V%, for the 
breathing mode, at corresponding frequency u >±. For fur¬ 
ther detail on the model see Ref. [201 . 

The inter-mode interaction is Hj = 2Kx a Xb ~ 
K(b A tf B + V A b B ), where K = Gm/oj m d 3 is the coupling 
rate with G the universal gravitational constant and d 
the distance between the mechanical resonators. One 
then diagonalizes the total Hamiltonian by applying the 
unitary operator U = exp{( 7 r/ 4 )( 6 y 4 &^ — S^S#)}. The 
effect of Gravitational coupling will appear in the free 
evolution part of Eq. ( |B5[ ) and after applying the above 
unitary transformation the equation reads 


Pge,eg — 


H ilUjtybj, Pge,eg\ 

3 =+,~ 

+ )Pge,eg 

.A e B 


+ i-^Pge,eg(b+ + b\ - b- - bl_) 

+ •£m i Pge,eg + £-qjPge,eg, 

J=+,- 1 = 1,2 


(Cl) 


where b± = (b A ± b B )/V2 are the normal annihilation 
operators and p ge , eg = U p geteg W. 

The initial conditions (the coupling signs to the qubits 
and the initial positions) are fixed such that only one nor¬ 
mal mode is excited and put into a superposition via her¬ 
alded technique. In realistic situation the normal mode 
splitting is expected to be small, and thus is approxi¬ 
mated to be A = w_ — w + ss 2Gm/ojd 3 [20] . In the first 
scenario we imagine that only the center-of-mass mode 
b + with frequency is excited by the proper constant 
coupling strength t o th e qubit: A e ^ A = — X e ,B = A e . This 
corresponds to Eq.(|B8|) and gives 


Tr„ 


Ap 


ge,eg(t)} = ±*e 27<,T Tr m+ {V + [\f2a e {r))p{t 0 )}. 


Here, V + (y/2a e ) = = V A {a e )V B {a e ). 

While for the second scenario where only breathing mode 
6 _ is excited with the proper coupling strength to the 
qubits: \ e ,A = A e ,B = A e . This corresponds to Eq.(B7l. 
Therefore, we have 


Tr m A ,m B {Pge,e g (t)} = ±\e 1 2 ^ T Tr m _{2?_[v / 2o! e (r)]p(to)}- 


oscillating with normal mode frequency . 


Therefore, in this representation we can think of a sin¬ 
gle massive system being in a harmonic potential oscillat¬ 
ing with frequency u >±. Environmental noise as described 
by (B3) is coupled locally to each resonator’s coordi¬ 


nate, and therefore qubit pair entanglement, according 
to the definition (B6), undergoes the same decoherence 


rate regardless of the collective mode. In the presence of 
large enough normal mode spitting , due to gravitational 
force between the mechanical resonators, the collective 
mechanical modes undergo different gravitational deco¬ 
herence rates. The difference in the decoherence gen¬ 
erates a gap in the amount of entanglement revival per 
single oscillation period. The effect of normal mode fre¬ 
quency splitting, i.e. the gap in the entanglement revival, 
becomes more pronounced for larger A e /w m . The deco¬ 
herence channel opens up by the gravitational interaction 
may captured by 


£gra vP = ^ r grav ) [®±, p]\ ■ (C2) 


This form of master equation leads to the suppression of 
superposition in the position coordinate of a collective 
mode. This effect can also be easily seen by rewriting 
the double commutator in the position space, 


-^gravt^i, [*£±,P(f)]] t r' graV ( a '± X±) P(x± > *£± ) ^)> 


describing in particular the spatial collapse of the oscil¬ 
lating mode due to gravitational decoherence. The grav¬ 
itational decoherence is proportional to the inverse of the 
normal mode frequency, i.e., Tg rav a l/cj^ [2tJ|. There¬ 
fore, unlike environmental decoherence gravitational de¬ 
coherence attribute different decoherence rates for the 
two collective modes which can be manifested in the gap 
in the amount of the entanglement revival between the 
two scenarios. This is a preliminary illustration of how 
non-Gaussian entangled state might be useful for probing 
quantum dynamical phenomena which are of purely grav¬ 
itational effect. This leads us to ask, can non-Gaussian 
mechanical entanglement provide us a sensitive detection 
of the genuine effect of gravitational decoherence? To es¬ 
timate this effect in real experimental situation one need 
to take into account all the practical limitations which is 
a challenging issue but not impossible. 
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